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Final  Report 

Luca  Dal  Negro 

Departmen  t  of  Electrical  and  Computer  Engineering  &  Photonics  Center, 

Boston  University,  Boston,  MA,  02215 

Personnel  in  place:  one  graduate  student  worked  on  this  project  since  the  beginning  and 
a  post-doctoral  associate.  The  graduate  student  (Jacob  Trevino)  will  graduate  in  the 
Spring  2013  based  on  the  results  of  this  research. 

Summary  of  results:  during  this  project  we  have  developed  novel  approaches  to  enhance 
light-matter  coupling  in  optical  devices  by  engineering  photonic -plasmonic  resonances  in 
aperiodically  ordered  nanostructures.  In  particular,  we  have  designed,  fabricated  and 
characterized  a  large  number  of  active  photonic  nanomaterials  and  structures  and 
demonstrated  unique  optical  properties  such  as  broadband  enhanced  local  field  intensity, 
scattering,  radiative  (i.e.,  light  emission)  and  nonlinear  responses  (i.e.,  second  harmonic 
generation,  nonlinear  refractive  index)  using  Si  compatible  materials  and  processing. 
Fabricated  metal-dielectric  devices  have  been  integrated  with  planar  Si  chips  for 
applications  of  aperiodic  order  to  light  emission,  optical  sensing,  on-chip  nonlinear 
generation,  solar  energy  conversion,  and  singular  optics.  The  fabrication  and 
experimental  characterization  of  materials  and  device  demonstrators  have  been  partnered 
with  the  development  of  efficient  computational  design  tools  based  on  semi-analytical 
multiple  scattering  theory,  which  allowed  the  rigorous  study  and  optimization  of  large- 
scale  aperiodic  media  for  the  first  time. 

Productivity  metrics:  The  many  results  generated  during  this  project  have  led  to  the 
publication  of  32  peer-reviewed  journal  publications,  4  book  chapters,  and  1  graduate 
level  textbook  (the  first  on  the  topic)  titled:  Optics  of  Aperiodic  Media,  which  will  soon 
be  published  by  Cambridge  University  Press  (see  the  complete  publication  list  below). 
The  vast  range  of  activities  sponsored  by  this  project  has  also  resulted  in  27  invited  talks, 
all  delivered  by  Dal  Negro  at  major  international  conferences  in  photonics  and  optical 
materials,  as  well  as  17  contributed  talks  delivered  by  students  and  postdocs  and  18 
invited  seminars  at  Universities  and  leading  Scientific  Institutions. 

Main  scientific  achievements:  design  of  aperiodic  structures 

We  have  designed  and  engineered  aperiodic  metal  nanoparticle  arrays  with  controllable 
Fourier  spectral  properties  that  interpolate  in  a  tunable  fashion  between  disordered 
random  systems  and  regular  periodic  structures. 

We  refer  to  this  general  class  of  artificial  metal-dielectric  materials  as  Deterministic 
Aperiodic  Nano  Structures  (DANS)[21].  We  have  demonstrated  that  these  novel 
aperiodic  media,  conceived  by  designing  spatial  frequencies  in  aperiodic  Fourier  space, 
give  rise  to  characteristic  scattering  resonances  and  localized  mode  patterns  that  enhance 
the  intensity  of  optical  fields  over  planar  surfaces  and  broad  frequency  spectra.  Moreover, 
our  work  on  aperiodic  plasmon  arrays  has  unveiled  for  the  first  time  the  distinctive 
interplay  between  photonic  diffraction  and  near  field  plasmonic  localization,  providing 
novel  opportunities  to  manipulate  light-matter  interactions  on  the  nanoscale. 
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By  engineering  photonic -plasmonic  coupling  at  multiple  wavelengths  in  aperiodic 
nanoparticle  arrays  we  have  experimentally  demonstrated  broadband  plasmon-enhanced 
scattering  across  the  entire  solar  spectrum.  By  utilizing  rigorous  electromagnetic  design 
tools  [1-3,11,19,23]  (i.e.,  multi-particle  generalized  Mie  theory,  Null  Field  method, 
Coupled  Dipole  Theory)  we  have  determined  the  relevant  engineering  design  rules  for  the 
control  and  enhancement  of  optical  fields  in  complex  aperiodic  media,  and  supported  our 
findings  by  genetic  and  particle  swarm  optimization  techniques  [3,19].  A  number  of  high 
impact  papers  have  been  published  on  these  topics  (see  complete  list  below).  Our  work 
has  shown  that  the  ability  to  engineer  aperiodic  structures  in  between  random  and 
periodic  media  is  important  to  engineer  intense  hot  spots  distributed  over  large  cross 
sectional  areas  of  the  arrays.  This  point  was  illustrated  by  Forestiere  et  al.  [1]  who 
systematically  investigated  the  near-field  spectra  and  far-field  scattering  response  of  Ag 
nanoparticle  arrays  generated  according  to  prime  numbers  distributions  in  two  spatial 
dimensions.  Using  rigorous  coupled-dipole  analysis  for  dipolar  nanoparticles,  this  study 
demonstrated  that  significantly  increased  local  field  intensity  over  broad  frequency 
spectra  can  be  obtained  by  designing  closely  packed  arrays  (i.e.,  large  particle  filling 
fractions)  and  with  a  large  density  of  spatial  frequencies,  captured  by  the  arrays  spectral 
flatness1  (SF). 
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Fig.  1.  (a)  Maximum  field  enhancement  versus  the  wavelength  for  an  isolated  Ag  nanosphere  (50nm 
radius)  and  for  periodic,  Fibonacci,  Thue-Morse,  and  Rudin-Shapiro  aperiodic  arrays  of  nano-spheres  with 
50  nm  minimum  interparticle  separation.  The  arrays  are  excited  by  a  circularly  polarized  plane  wave  at 
normal  incidence,  (b)  Values  of  maximum  field  enhancement  versus  the  spectral  flatness  (SF)  x  filling 
fraction  (FF)  product  for  the  different  arrays  indicated  in  the  figure. 

This  point  is  also  illustrated  in  Fig.  1(a),  where  we  show  the  calculated  near- field  spectra 
for  Ag  nanosphere  arrays  with  periodic,  Fibonacci,  Thue-Morse,  and  Rudin-Shapiro 
sequences.  The  case  of  a  single  Ag  nanosphere  with  50nm  radius  is  also  shown  for 
comparison.  To  better  describe  the  influence  of  both  the  polarization  states  of  the  incident 
field,  the  arrays  were  excited  by  a  circularly  polarized  plane  wave  at  normal  incidence. 
Figure  1(b)  demonstrates  the  scaling  behavior  of  the  maximum  hot  spot  field  intensity, 
probed  in  the  plane  of  the  arrays,  for  a  large  number  of  aperiodic  deterministic  structures 


!The  spectral  flatness  (SF)  is  a  digital  signal  processing  parameter  that  measures  how  spectrally  diffused  a 
signal  is.  In  the  case  of  plasmonic  structures,  the  arrays  are  considered  as  digitized  2D  spatial  signals  and  the 
SF  is  calculated  by  dividing  the  geometric  mean  and  the  arithmetic  mean  of  their  Fourier  power  spectra. 
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(named  in  the  figure)  as  a  function  of  the  product  of  the  arrays  filling  fraction  (FF)  and 
their  Fourier  spectral  flatness  (SF).  This  plot  demonstrates  that  only  aperiodic  structures 
which  possess  both  large  spectral  flatness  SF  and  nanoparticle  filling  fraction  result  in 
strong  plasmonic  scattering  and  near-field  plasmonic  localization.  The  rationale  for  this  is 
that  plasmon  waves  couple  very  strongly  in  the  near-field  regime  at  very  short  distances, 
thus  motivating  the  need  of  high  nanoparticles  packing  fractions.  However,  for  a  given 
particle  density,  aperiodic  arrays  featuring  a  large  number  of  spatial  frequencies  can 
match  in-plane  scattering  processes  resulting  in  efficient  multiple  scattering  in  the  array 
plane,  and  boosting  the  field  enhancement  even  further.  The  coexistence  of  different 
electromagnetic  coupling  regimes  at  multiple  length  scales  is  at  the  origin  of  the  superior 
field  enhancement  and  localization  observed  in  several  aperiodic  plasmonic  structures. 
Due  to  the  increased  structural  disorder  (spectral  flatness),  the  nanoparticles  in  the 
aperiodic  arrays  are  strongly  coupled  in  both  the  plasmonic  near  field  regime  and  the 
photonic  diffractive  one  (long-range  coupling),  resulting  in  strong  in-plane  multiple  light 
scattering.  On  the  other  hand,  in  the  case  of  periodic  array  the  nanoparticles  are 
prevalently  coupled  in  the  near-field  regime,  with  strongly  reduced  phase  modulation 
across  the  array  plane.  In  this  case,  extended  plasmonic  structural  resonances,  analogous 
to  the  Fabry-Perot-type  modes  in  finite-size  structures,  are  formed  and  decrease  the 
overall  field  intensity  enhancement. 
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Fig.  2.  The  color-maps  show  the  cumulative  distributions  of  field  enhancement  (CDFE)  (logarithmic  scale) 
versus  wavelength  (x-axis)  and  field -enhancement  (y-axis)  for  (a)  Periodic,  (b)  Fibonacci,  (c)  Thue-Morse, 
(d)  Rudin-Shapiro.  The  arrays  consist  of  50  nm  radius  Ag  spheres  with  minimum  edge-to-edge  separations 
of  25  nm.  They  are  excited  by  a  circularly  polarized  plane  wave  at  normal  incidence. 

It  is  finally  to  be  noticed  that  aperiodic  arrays  have  been  found  to  perform  much  better 
than  closely  packed  periodic  ones  despite  their  reduced  values  of  filling  fraction. 
Therefore,  the  results  in  Fig.  1  demonstrate  clearly  that  closely-packed  aperiodic 
plasmonic  arrays  with  a  large  density  of  spatial  frequency  are  needed  in  order  to  enhance 
the  hot  spots  intensity  over  a  broader  frequency  range  compared  to  optimized  periodic 
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and  quasiperiodic  structures.  Another  important  aspect  of  aperiodic  plasmonic  arrays  that 
we  have  discovered  is  related  to  the  fraction  of  the  total  surface  area  covered  by  strong 
plasmonic  fields.  In  plasmonic  sensing  technology,  the  control  of  the  areal  density  of 
enhanced  fields  on  planar  chips  is  of  fundamental  importance.  In  order  to  quantitatively 
understand  this  aspect,  we  have  studied  the  fraction  of  the  total  area  of  the  arrays  covered 
by  plasmon  enhanced  fields  with  values  greater  then  a  fixed  threshold  value. 
Mathematically,  this  can  be  quantified  by  introducing  the  cumulative  distribution  of  field 
enhancement  (CDFE),  discussed  in  Ref.[l] 

In  Fig.  2  we  show  the  calculated  CDFE  for  (a)  periodic,  (b)  Fibonacci,  (c)  Thue-Morse, 
and  (d)  Rudin-Shapiro  two-dimensional  arrays,  respectively.  The  CDFE  function 
describes  the  fraction  of  the  total  area  of  these  arrays  that  is  covered  by  enhanced  fields 
with  values  greater  than  the  value  specified  in  the  vertical  axis.  Although  the  CDFE  does 
not  provide  any  information  about  the  size  of  hot  spots,  it  gives  a  quantitative  measure  of 
the  spatial  distribution  of  enhanced  fields,  at  each  wavelength,  with  respect  to  the  total 
surface  of  the  array.  The  results  in  Fig.  2  demonstrate  that  aperiodic  arrays  with  large 
spectral  flatness  and  particle  filling  fraction  support  enhanced  field  states  that  are 
spatially  distributed  over  larger  array  areas  compared  to  periodic  plasmonic  structures, 
which  is  a  very  important  attribute  for  the  engineering  of  scattering-based  sensors  (e.g., 
SERS  substrates)  and  planar  plasmonic  devices.  However,  we  should  notice  here  that 
aperiodic  designs  with  dense  Fourier  spectra  come  at  the  additional  cost  of  a  larger 
system’s  size  compared  to  narrow-band  periodic  or  multi -periodic  structures,  ultimately 
requiring  engineering  trade-offs  between  the  intensity  enhancement,  the  resonant 
frequency  bandwidth,  and  the  total  size  of  plasmonic  devices. 

The  characteristic  behavior  of  aperiodic  nanoplasmonic  structures  summarized  in  Figs.  1 
and  2  follows  from  the  large  number  of  supported  photonic  modes  available.  The  large 
spectral  density  of  photonic  resonances  distinctive  of  aperiodic  media  couple  to  sub¬ 
wavelength  localized  plasmon  modes  over  a  broad  frequency  range,  controlled  by  the 
structural  complexity  of  the  array.  Broadband  hot  spots  intensity  enhancement  with 
aperiodic  plasmonic  structures  is  therefore  made  possible  by  controlling  long-range 
electromagnetic  coupling  in  multi-scale  arrays  of  resonant  nanoparticles  with  positioned 
fluctuations. 

Another  important  featured  that  we  discovered  in  multi-scale  electromagnetically  coupled 
aperiodic  arrays  of  metallic  nanoparticles  is  the  distinctive  scaling  of  their  hot  spots 
intensity  with  respect  to  the  system’s  size,  or  the  total  number  of  particles  in  the  array. 
This  is  best  illustrated  by  the  analytical  multiple  scattering  results  shown  in  in  Fig.  3  for 
Au  spheres.  Additionally  we  illustrate  the  field  intensity  distribution  in  the  plane  of  the 
arrays  for  periodic  square  arrays,  Fibonacci,  Thue-Morse,  and  Rudin-Shapiro  structures. 
All  the  structures  are  illuminated  at  normal  incidence  by  a  plane  wave  at  785nm.  The 
distinctive  interplay  between  near-field  plasmonic  coupling  and  long-range  multiple 
scattering  in  aperiodic  structures  is  clearly  displayed  in  Figs.  3(b-d)  where  we  show  the 
electric  field  spatial  distribution  at  the  wavelengths  of  the  enhancement  peak.  The  highly 
inhomogeneous  field  distributions  in  Figs.  3(b-d)  demonstrate  the  importance  of  long- 
range  multiple  scattering  effects  in  the  plane  of  the  array,  which  couple  all  the  particle 
clusters  in  the  aperiodic  arrays.  Moreover,  when  scaling  up  the  size  of  the  arrays  by 
increasing  the  particles  number  (Figs.  3e),  new  configurations  of  local  particle  clusters  in 
aperiodic  arrays  appear  separated  by  wavelength- scale  distance,  thus  increasing  the  total 
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number  of  spatial  frequencies  in  the  plane  and  enhancing  the  maximum  hot  spots 
intensity  in  a  size-dependent  fashion. 


Fig.  3.  Generalized  Mie  Theory  (GMT)  calculations  of  electromagnetic  field  scattered  by  plasmonic  arrays 
of  spherical  Au  nanoparticles  arranged  according  to  (a)  a  periodic  (b)  Fibonacci  (c)  Thue-Morse  (d)  Rudin- 
Shapiro  array.  All  the  particles  in  the  arrays  have  a  radius  of  lOOnm  and  a  minimum  separation  of  25nm. 
(e)  Scaling  of  the  calculated  intensity  enhancement  for  the  different  arrays  as  a  function  of  the  number  of 
particles. 

The  size  dependent  nature  of  the  optical  response  of  aperiodic  systems  is  a  direct 
manifestation  of  multiple  scattering  in  the  mesoscopic  regime.  On  the  other  hand,  no 
photonic  coupling  occurs  in  subwavelength  packed  periodic  array,  as  evidenced  by  the 
strongly  reduced  phase  modulation  across  the  array  plane  in  Fig.  3(a),  and  more 
delocalized  plasmonic  modes  are  formed  across  the  entire  periodic  structure  with  reduced 
hot  spot  intensity.  The  lack  of  photonic-type  coupling  in  closely  packed  periodic  arrays 
prevents  the  onset  of  size-dependent  photonic -plasmonic  resonances  and  the  resulting 
field  enhancement  effects,  making  the  maximum  hot  spot  intensity  almost  insensitive  to 
the  overall  periodic  array  size,  as  shown  in  Fig.  3(e). 

We  conclude  by  noticing  that  the  mesoscopic  regime  of  aperiodic  plasmon  arrays  has 
direct  implications  for  optical  device  engineering.  In  fact,  in  addition  to  the  particles 
composition/morphology  and  array  geometry,  we  have  demonstrated  that  the  size  of 
aperiodic  arrays  can  be  tailored  in  order  to  largely  enhance  the  near-field  hot  spots 
intensity. 

Nanofabrication  of  aperiodic  nanostructures 

In  contrast  to  random  media,  DANS  can  be  specifically  tailored  and  fabricated  using 
conventional  nanolithographic  techniques  such  as  electron  beam  lithography  (EBL)  or 
Focused  Ion  Beam  (FIB)  milling  followed  by  standard  metal  deposition  and  etching  steps. 
Within  this  project,  our  group  has  recently  developed  a  flexible  process  flow,  sketched 
below,  for  the  nanofabrication  of  arbitrary  arrays  of  metal  nanoparticles  of  interest  to 
nanoplasmonic  applications.  Moreover,  high-throughput  fabrication  was  demonstrated 
recently  by  Lin  et  al.  [29]  who  combined  this  approach  with  a  reusable  transfer  imprint 
technique.  Aperiodic  nanoparticle  arrays  based  on  noble  metals,  typically  Au  and  Ag,  are 


Deterministic  Aperiodic  Structures  for  on-chip  nanophotonics 
and  nanoplasmonics  device  applications,  Award  FA9550-10-1-0019 


typically  fabricated  on  quartz  substrates  using  a  10  nm  layer  Indium  Tin  Oxide  (ITO)  to 
provide  conduction. 
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Fig.  4.  Schematics  of  the  fabrication  process  flow  for  the  generation  of  aperiodic  plasmonic  nanoparticle 
arrays  (left)  and  nano-hole  arrays  (right).  The  SEM  pictures  show  a  Rudin-Shapiro  arrays  of  Au 
nanoparticles  (left)  and  a  Gaussian  prime  array  of  nano-holes  in  a  quartz  substrate.  The  particle/hole  radius 
is  lOOnm  and  the  minimum  interparticle  separation  is  200nm. 


A  1 80-nm-thick  layer  of  PMMA  (PolyMethylMethAcrylate)  is  then  spin  coated  on  top  of 
the  cleaned  substrate.  Subsequently,  the  DANS  patterns  are  defined  using  a  Zeiss  SUPRA 
40VP  SEM  equipped  with  a  Raith  Beam  Blanker  and  NPGS  for  nanopatterning.  After 
developing  the  resist  in  a  1:3  solution  of  MIBK  (Methyl  IsoButyle  Ketone)  and  IPA 
(Isopropanol),  a  ~  30nm  thick  Au/Ag  film  is  deposited  on  the  patterned  surface  by 
electron-beam  evaporation.  Finally,  a  liftoff  process  is  performed  using  acetone,  resulting 
in  the  definition  of  the  targeted  metal  nanoparticle  arrays. 

Within  the  same  general  process  flow,  nano-perforated  metal/dielectric  films  can  also  be 
obtained  using  by  a  Reactive  Ion  Etching  (RIE)  step  immediately  after  the  EBL  writing. 
The  process  flow  for  the  fabrication  of  both  metallic  nanoparticle  arrays  and  nano-hole 
patterns  is  illustrated  in  Fig.  4.  Typical  dimensions  of  each  fabricated  array  are  about 
100pm  x  100pm  for  nanoparticles  with  a  diameter  of  200nm,  30nm  tall  and  variable 
separations  that  can  range  in  between  25nm  and  400nm,  depending  on  the  specific  DANS 
geometry,  device  applications,  and  nanolithographic  setup. 

As  an  example  of  fabricated  DANS,  we  show  in  Figure  5  the  Scanning  Electron 
Microscopy  (SEM)  pictures  of  arrays  of  Au  nanoparticles  arranged  in  Fibonacci  (a), 
Thue-Morse  (b),  Rudin-Shapiro  (c),  and  co-prime  (d)  pattern  geometries.  The  Au 
particles  are  cylindrical  in  shape  and  their  height,  as  characterized  by  Atomic  Force 
Microscopy  (AFM)  and  SEM,  was  found  to  be  h=30nm.  All  the  particles  have  a  circular 
diameter  of  d=200nm  and  a  minimum  interparticle  separation  a=25nm.  Notice  however, 
that  additional  length  scales  are  present  in  aperiodic  structures,  extending  to  dimensions 
comparable  to  the  wavelength  of  light  in  the  optical  regime. 
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Fig.  5.  SEM  pictures  of  (a)  Fibonacci,  (b)  Thue-Morse,  (c)  Rudin-Shapiro,  (d)  Coprime  Au  nanoparticle 
array.  The  individual  particle  sizes  are  150nm,  and  the  minimum  interparticle  separations  in  the  arrays 
shown  are  25nm. 

To  demonstrate  the  flexibility  of  the  DANS  fabrication  approach  described  above,  we 
show  in  Fig.  6  additional  examples  of  Au  nanoparticle  arrays  fabricated  with  various 
types  of  deterministic  aperiodic  order  on  Si  substrates.  These  structures  are  a  Danzer  (a), 
Pinwheel  (b),  body-centered  Pinwheel  (c)  and  the  three  most  investigated  types  of  Vogel 
spiral  arrays,  namely  the  ai-spiral  (d),  the  golden  angle  spiral  (e)  and  the  a2-spiral  (f).  We 
can  clearly  appreciate  from  Figures  5  and  6  the  quality  of  the  nanofabricated  arrays, 
which  excellently  match  the  designed  geometrical  patterns  over  large  device  areas. 
However,  current  nanoscale  writing  techniques,  such  as  EBL,  focused  ion  beam 
lithography  (FIB),  and  scanning  probe  microscopy  (SPM),  suffer  from  high  operating 
costs  and  low  throughput.  In  order  to  efficiently  scale  the  dimensions  of  aperiodic 
plasmonic  arrays  to  larger  sizes  in  a  cost-effective  manner,  novel  fabrication  methods  are 
currently  investigated.  We  notice  that  shadow-mask  patterning  techniques,  such  as 
stencil-based  methods,  allow  for  the  large-scale  fabrication  of  plasmonic  nanostructures 
but  inherently  suffer  from  edge  blurring  and  cannot  produce  plasmonic  films  perforated 
with  nano-holes,  which  play  an  important  role  in  nanoplasmonics. 
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Fig.  6.  SEM  pictures  of  (a)  Danzer,  (b)  Peenwheel,  (c)  Body  centered  pinwheel,  (d)  ai  Spiral  (137.3°),  (e) 
golden  Angle  Spiral,  (f)  a2  Spiral  (137.6°)  Au  nanoparticle  array.  The  individual  particle  sizes  are  200nm, 
and  the  minimum  interparticle  separations  in  the  arrays  shown  are  50nm. 

To  solve  these  problems,  we  have  demonstrated  a  scalable  and  cost-effective  direct 
transfer  nanofabrication  technique  that  utilizes  a  hard  mold  master  and  an  inexpensive, 
commercially  available  flip-chip  bonder,  for  the  fabrication  of  large-scale  metallic 
nanoparticles  on  polymer  substrates  and  perforated  metallic  membranes  (nano-hole 
arrays)  atop  silk  fibroin  films  [29]. 

The  process  flow  for  the  transfer  imprint  of  plasmonic  nano-dots  and  nano-holes,  begins 
with  the  fabrication  of  the  reusable  master  molds.  We  illustrate  in  Fig.  7  the  successful 
printing  of  nano-dots  on  silk  fibroins  substrates  and  nano-hole  arrays  on  a  metal  film.  In 
the  case  of  the  nano-dot  transfer  process  (Fig.  7a-d),  the  desired  geometry  is  fabricated 
into  a  Si  mold  consisting  of  nanopillar  arrays,  while  the  nano -hole  process  requires  a 
mold  containing  nano-holes,  as  shown  in  Fig.  7  (e-h).  The  fabrication  of  the  nano-hole 
master  proceeds  via  EBL  writing  with  a  260nm-deep  RIE  step,  using  the  PMMA  as  an 
etch  mask.  The  remaining  PMMA  is  removed  by  hot  acetone  bath,  resulting  in  the  Si 
nano-hole  master.  The  Si  master  is  first  treated  with  a  silanizing  agent  to  reduce  the 
adhesion  of  the  Au  to  the  Si  surface.  This  surface  treatment  enables  a  higher  yield  in 
pattern  transfer  of  the  Au  to  the  silk  film  in  the  subsequent  steps.  The  process  flow 
continues  with  the  deposition  of  a  35nm-thick  e-beam  evaporated  gold  (Au)  film,  as 
shown  in  Fig.  7  (f). 
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Fig.  7.  Process  flow  of  transfer  nanoimprint  of  plasmonic  nano-dots  (a-d)  and  nano-hole  arrays  (e-h)  using 
reusable  masters,  (i)  Transferred  plasmonic  at  spiral  nanodot  array  on  silk  film  from  Si  master.  (1) 
Transferred  plasmonic  ai  spiral  nano-hole  array  on  silk  film  from  Si  master. 

The  Au  coated  master  is  now  ready  for  transfer  imprinting  the  plasmon  membrane  on  the 
silk  substrate  (Fig.  7g-h).  A  commercial  flip-chip  bonder  (Smart  Equipment  Technology 
FC150)  was  employed  to  transfer  imprint  from  the  fabricated  master  mold.  The  flip-chip 
bonder  is  used  to  align  and  bond  one  or  more  chips  onto  a  substrate  using  pressure  and 
heat.  The  transfer  imprint  process  begins  by  heating  the  Au  covered  master  to  90°C,  when 

9 

it  is  then  pressed  onto  the  surface  of  the  silk  layer  with  a  force  of  50kg  (area  of  1 6 mm") 
for  5min.  During  the  imprinting,  the  plasmonic  nanostructures  bind  to  the  surface  of  silk 
layer.  Upon  completion  of  the  bonding  cycle,  the  master  is  removed,  leaving  the  Au 
embedded  on  the  surface  of  the  silk  layer  as  shown  in  Fig.  7(h)  for  the  case  on  nano- 
perforated  membranes.  Nano -imprinted  plasmonic  hole  arrays  with  a  representative 
Vogel  spiral  geometry  [8]  are  shown  as  an  example  in  Fig.  7  (1).  A  similar  process  has 
been  successfully  developed  to  imprint  arrays  of  metallic  nano-dots  atop  polymer  films. 
In  this  case,  a  Si  mold  consisting  of  nanopillar  arrays  needs  to  be  utilized,  as  illustrated  in 
Figs.  7(a-d).  The  nanodot  master  fabrication  proceeds  with  the  deposition  of  35nm-thick 
chrome  (Cr)  layer  via  electron  beam  (e-beam)  evaporation,  followed  by  a  lift-off  process 
in  heated  acetone.  The  resulting  Cr  nanocylinders  serve  as  a  etch  mask  for  a  230nm-deep 
anisotropic  reactive  ion  etch  (RIE),  forming  the  desired  nanopillar  array.  The  residual  Cr 
layer  is  then  removed  by  wet  chemical  etch.  More  details  on  these  novel  high  throughput 
and  scalable  fabrication  processes  can  be  found  in  Ref  [29]. 

Reusable  transfer  imprint  enables  large-scale  replication  of  arbitrarily  complex 
nanostructures  with  deep  sub-wavelength  details  down  to  30nm  with  a  high  throughput. 
The  inexpensive  scalability  of  aperiodic  nanoplasmonic  structures  over  large  areas  is  of 
great  importance  for  device  engineering  as  it  offers  the  opportunity  to  greatly  reduce 
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fabrication  costs  and  to  develop  aperiodic  substrates  into  a  mature  technology  for  cost- 
effective  device  applications. 

Light  emission  enhancement 

We  designed  and  fabricated  Au  nanoparticle  arrays  with  quasi-periodic  (Fibonacci)  order 
and  coupled  them  to  light  emitting  Er-doped  silicon  nitride  materials  in  order  to 
investigate  plasmon  enhancement  of  1.55  microns  radiation  [2]. 

A  3.6  times  enhancement  of  the  photoluminescence  (PL)  intensity  accompanied  by  a 
reduction  of  the  Er3+  emission  lifetime  at  1.54pm  (see  Fig.  8)  has  been  observed  and 
explained  by  the  coupling  of  light  emission  with  the  broadband  near-infrared  structural 
resonances  supported  by  the  Fibonacci  structure.  Photonic -plasmonic  hybrid  Fibonacci 
modes  were  identified  in  a  comparative  transmission  experiment  between  periodic 
plasmonic  gratings  and  Fibonacci  quasi-periodic  structures  by  the  spectral  positions  of 
transmission  minima,  in  a  large  range  of  fabricated  structures.  The  scaling  behavior  of  the 
emission  enhancement  in  Fibonacci  and  periodic  gratings  was  discussed  in  relation  to  the 
modifications  of  the  photonic  Local  Density  of  States  (LDOS)  at  the  Er  emission 
wavelength.  The  strongly  reduced  frequency  sensitivity  of  the  photonic-plasmonic 
scattering  resonances  in  Fibonacci  arrays  prevented  the  detuning  from  the  Er  emission 
wavelength,  and  increased  Er  emission  across  its  entire  spectrum. 
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Fig.  8.  Demonstration  of  light  emission  enhancement  from  Erbium  atoms  coupled  to  deterministic 
aperiodic  plasmonic  arrays  of  Au  nanoparticles  (200nm  diameter)  (a)  SEM  picture  of  quasiperiodic 
Fibonacci  Au  nanoparticle  array  fabricated  atop  a  light  emitting  Er:SiNx  substrates  of  80nm  thickness,  (b) 
Integrated  PL  enhancement  in  periodic  (black)  and  Fibonacci  (red)  arrays  of  various  interparticle 
separations.  The  inset  shows  representative  PL  spectra  of  periodic  (black),  Fibonacci  (red)  arrays  and  of  the 
unpatterned  (green)  region,  (c)  PL  decay  time  of  Er  atoms  through  unpatterned  substrate  (black)  and 
Fibonacci  arrays  with  varying  interparticle  separations  indicated  in  the  legend,  (d)  PL  spectra  excited  at 
488nm  of  periodic  and  aperiodic  nanoparticle  arrays  (50nm  min  interparticle  separation). 
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Moreover,  when  decreasing  the  interparticle  separation  in  the  Fibonacci  arrays,  the 
emission  intensity  was  found  to  increase  along  with  the  emission  rate,  while  an  opposite 
behavior,  indicative  of  non-radiative  losses,  was  observed  in  periodic  grating  structures. 
The  coupling  of  light  emission  to  the  broadband  scattering  resonances  of  quasi-periodic 
Fibonacci  arrays  makes  these  systems  particularly  attractive  for  radiative  rate  engineering 
applications  that  require  light  extraction/enhancement  over  broad  frequency  spectra. 

We  also  developed  an  alternative  approach  to  light  emission  enhancement  based  on  the 
engineering  of  circular  light  scattering  in  plasmonic  Vogel’s  spirals  [8].  We  demonstrated 
that  this  can  result  in  broadband  light  emission  enhancement  in  active  thin  films  coupled 
to  deterministic  aperiodic  structures.  Within  this  project  we  have  investigated  the  impact 
of  the  inhomogeneous  distribution  of  local  spatial  frequencies  of  Vogel  spirals  on  their 
emission  patters.  In  particular,  we  prepared  a  dye  polymer  solution  by  dissolving 
common  laser  dye  molecules  of  DCM  (Exciton  Inc.),  in  toluene. 


A  /=\ 

CH—Z  >-N(CH.), 


Laser  Die:  DCM  (Exciton,  Inc.) 


lOOnm  t 


Dye  Polymer 
■  Spiral  Array 


Quartz  Substrate 


550  600  650  700  750  800  850 


Wavelength  (nm) 


I  Au  nanoparticles 
2  Dye  polymer 
|  Quartz  substrate 


Fig.  9.  (top)  Schematic  cross-section  of  Au  spiral  nanoparticle  array  on  a  quartz  substrate,  coated  with 
DCM  doped  PMMA.  Nanoparticles  are  cylindrical  in  shape  with  200nm  diameter  and  30nm  in  thickness 
with  array  diameter  of  100pm.  (a-d)  CCD  images  of  light  emission  from  a  DCM  dye  layer  (lOOnm  thick) 
deposited  onto:  (a)  homogeneous  quartz  substrate  (b)  dye  emission  from  scattered  light  of  a2-spiral  (27,778 
particles)  with  laser  positioned  at  the  center  of  the  spiral,  (c)  and  (d)  dye  emission  from  scattered  light  of 
a2-spiral  with  laser  positioned  off  center  of  the  spiral,  (e)  Photoluminescence  spectra  of  doped  polymer 
form  samples  coupled  to  a2-spiral  arrays  with  varying  average  interparticle  separations  (310nm,  570nm, 
900nm  and  un-coupled  film  for  reference). 


The  dilute  solution  was  then  mixed  with  Polymethylmethacrylate  (PMMA),  spun  onto 
samples  and  cured,  resulting  in  lOOnm  thick  films  of  laser  dye  doped  PMMA.  This 
particular  laser  dye  has  maximum  absorption  at  480nm  and  an  emission  peak  at  640nm, 
which  overlaps  the  scattering  resonances  of  the  investigated  ot2-spiral,  leading  to 
broadband  emission  enhancement,  shown  in  Fig.  9.  Moreover,  we  demonstrate  the  ability 
to  dramatically  modify  the  angular  emission  of  the  sample  by  imaging  the  fluorescence  in 
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transmission  under  for  different  excitation  conditions  (Figs.  9a-d).  In  particular,  the 
sample  was  pumped  by  positioning  the  laser  spot  at  different  locations  onto  the  doped 
PMMA  substrate  at  normal  incidence  (focused  through  a  lOx  objective)  with  a  laser 
diode  at  480nm  and  the  emitted  light  was  collected  in  transmission  configuration  through 
the  substrate  using  a  lens  of  100mm  focal  length,  and  imaged  by  a  CCD  camera.  In  order 
to  capture  only  the  emission  patterns,  the  pump  laser  light  was  blocked  by  a  514nm  high- 
pass  filter.  An  identically  prepared  emitting  layer  was  also  coated  on  unpattemed  quartz 
for  reference.  Figure  9  (a)  shows  the  CCD  image  of  the  fluorescence  collected  in 
transmission  through  the  reference  sample  (i.e.,  with  no  spiral  pattern),  which  indicates 
that  the  fluorescence  is  spatially  confined  to  the  pumped  region  in  the  absence  of 
scattering  structures.  On  the  other  hand,  when  pumping  the  samples  with  the  plasmonic 
arrays,  the  fluorescence  spreads  laterally  in  the  plane  of  the  array  and,  when  the  sample  is 
symmetrically  pumped  through  its  center,  a  significant  fraction  of  the  fluorescence  is 
emitted  along  multiple  directions  due  to  the  isotropic  character  of  the  Fourier  space  (i.e., 
circular  light  scattering),  as  demonstrated  in  Figs.  9  (b).  Moreover,  when  the  position  of 
the  laser  pumping  spot  is  slightly  misplaced  from  the  center  of  the  sample  (in  the 
horizontal  direction)  by  approximately  25  pm,  the  angular  distribution  of  the  radiation 
changes  dramatically  due  to  the  very  inhomogeneous  distribution  of  local  spatial 
frequencies  associated  to  the  surface  of  aperiodic  spirals  with  circularly-symmetric 
Fourier  space,  as  shown  in  Figs.  9(c-d).  These  results  demonstrate  that  aperiodic  arrays 
with  Vogel  spiral  geometry  can  be  utilized  to  enhance  and  manipulate  light  emission  in 
planar  structures.  While  more  work  is  required  to  fully  leverage  aperiodic  spiral  order  for 
active  optical  devices,  we  believe  that  this  work  motivates  the  development  of  novel 
optical  devices  for  polarization  insensitive,  enhanced  light-matter  coupling  on  planar 
surfaces,  such  as  plasmonic  photodetectors  and  thin-film  solar  cells,  which  will  be  briefly 
reviewed  next. 

Thin-film  solar  and  photodetectors  enhancement 

In  order  to  broaden  the  spectral  region  of  light  absorption  enhancement,  it  is  crucial  to 
engineer  aperiodic  nanoparticle  arrays  with  a  higher  density  of  spatial  frequencies  and  a 
high  degree  of  rotational  symmetry  without  resorting  to  uncontrollable  random  systems, 
which  have  only  limited  engineering  appeal.  Recent  studies  have  proposed  to  utilize 
plasmonic  arrays  with  aperiodic  quasicrystal  structures,  such  as  Penrose  lattices,  which 
exhibit  non-crystallographic  rotational  symmetries.  Such  arrays,  by  virtue  of  their  higher 
degree  of  rotational  symmetry  compared  to  traditional  periodic  structures,  give  rise  to 
enhanced  scattering  along  multiple  directions  and  over  a  broader  wavelength  range, 
producing  coherent  wave  interactions  among  an  increased  number  of  neighboring 
particles.  There  is  currently  the  need  to  engineer  densely-packed  arrays  of 
nanoparticles/nano-holes  with  continuous  circular  symmetry  providing  polarization 
isotropy  and  efficient  planar  light  diffraction  within  targeted  spectral  ranges.  In  addition, 
such  effects  should  occur  regardless  of  the  direction  of  incident  radiation,  thus 
eliminating  the  need  to  costly  tracking  schemes  in  solar  cells  applications. 

In  this  project,  we  have  recently  proposed  and  utilized  plasmonic  aperiodic  Vogel  spiral 
arrays  as  a  viable  strategy  to  achieve  wide-angle  light  scattering  for  broadband  and 
polarization  insensitive  absorption  enhancement  in  thin-film  Si  solar  cells  [22].  Our 
group  has  previously  demonstrated  that  aperiodic  Vogel  spiral  arrays  feature  nearly 
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continuous  azimuthal  symmetry  in  Fourier  space  [8,18,25]  and,  when  normally 
illuminated,  satisfy  the  Rayleigh  condition  for  planar  light  scattering  over  broad  and 
controllable  frequency  bands,  irrespective  of  the  incident  polarization  of  light.  A  simple 
scalar  Fourier  optics  picture  already  suggests  that  polarization-insensitive  large-angle 
scattering  of  incident  radiation  occurs  in  GA  spiral  arrays  at  frequencies  matching  the 
radial  position  of  the  scattering  ring  in  reciprocal  space.  This  property  is  ideal  to 
manipulate  light  trapping  in  thin-film  Si  solar  cells.  Following  this  approach,  we  designed 
and  fabricated  GA  arrays  of  Au  nanoparticles  atop  ultra-thin  film  (i.e.,  50nm-thick 
absorbing  amorphous  Si  layer)  Silicon  On  Insulator  (SOI)  Schottky  photodetector 
structure,  as  sketched  in  Fig.  10,  and  we  demonstrated  experimentally  larger  photocurrent 
enhancement  in  the  600nm-950nm  spectral  range  compared  to  optimized  nanoparticle 
gratings.  The  relation  between  the  spatial  Fourier  spectrum  of  GA  arrays  of  Au 
nanoparticles  and  the  large  angular  distribution  of  scattered  radiation  in  the  forward 
scattering  hemisphere  has  been  rigorously  investigated  by  calculating  the  angular 
radiation  diagrams  using  the  Coupled  Dipole  Approximation  (CDA)  methods  for 
particles  with  ellipsoidal  shape. 


Fig.  10.  (top)  Schematics  of  the  device  cross-section  of  the  SOI  Schottky  photo -detector  with  plasmonic 
arrays  integrated  onto  the  absorbing  surface.  (a,b)  Calculated  radiation  diagrams  as  a  function  of  the 
inclination  angle  for  a  GA  spiral  (a)  and  a  periodic  array  (b)  for  three  different  wavelengths,  namely  480nm 
(Blue),  520nm  (Green),  and  650nm  (Red),  (c)  Spiral  array  photocurrent  with  average  center  to  center 
spacing  of  425nm  (red)  and  empty  neighbor  references  cells  (black),  (d)  Integrated  photocurrent 
enhancement  ratio  for  GA  spiral  (red)  and  periodic  (black  dashed)  arrays  of  different  center  to  center 
particle  spacing. 
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This  CDA  approach  is  particularly  suited  to  efficiently  treat  large-scale  plasmonic 
systems  of  small  and  well-separated  nanoparticles,  and  it  has  been  previously  validated 
against  semi-analytical  multiple  scattering  methods.  In  our  work,  all  nanoparticles  were 
modeled  by  oblate  spheroids  with  lOOnm  diameter  and  30nm  height.  Moreover,  the 
arrays  are  embedded  in  Si  and  normally  excited  by  a  linearly  polarized  plane  wave.  The 
parameter  of  interest  for  the  understanding  of  the  angular  scattering  properties  is  the 
differential  scattering  cross  section,  which  describes  the  angular  distribution  of 
electromagnetic  power  density  scattered  at  a  given  wavelength  within  a  unit  solid  angle 
centered  around  an  angular  direction  (0,cp)  per  unit  incident  irradiance.  In  the  case  of 
arrays  composed  of  dispersive  metal  nanoparticles,  the  power  scattered  from  a  particular 
structure  is  generally  a  function  of  both  the  geometrical  parameters  of  the  array  and  the 
wavelength  of  the  incident  radiation.  Full  information  on  angular  scattering  is  thus 
captured  by  calculating  the  averaged  differential  scattering  cross  section,  where  the 
average  is  performed  on  the  azimuthal  angle  (p  and  the  scattered  intensity  is  normalized  to 
the  maximum  value  (i.e.,  forward  scattering  peak).  By  plotting  the  azimuthally  averaged 
differential  scattering  cross  section  versus  the  inclination  angle,  we  obtain  the  radiation 
diagrams  of  the  arrays.  In  Figs.  10(a,b)  we  show  (plotted  in  dB  scale)  the  calculated 
radiation  diagrams  for  a  periodic  grating  structure  and  the  optimized  GA  arrays  at  three 
different  wavelengths  A|j=480nm,  /<;=550nm  and  Ar=6  1  Onm  (i.e  corresponding  to  the 
blue,  green,  and  red  colors),  respectively.  Differently  from  the  well-known  case  of 
periodic  structures,  where  the  scattered  radiation  is  preferentially  redistributed  along  the 
directions  of  coherent  Bragg  scattering  (Fig.  10a),  the  radiation  diagram  of  GA  arrays  is 
significantly  broadened  at  large  angles  (i.e.  >30°)  for  all  the  investigated  wavelengths, 
demonstrating  broadband  wide-angle  scattering  behavior  (Fig.  10b). 

By  combining  experimental  absorption  enhancement  and  photocurrent  measurements 
with  CDA  and  full-vector  3D  Finite  Difference  Time  Domain  (FDTD)  simulations,  we 
recently  demonstrated  that  broadband  wide-angle  scattering  in  GA  spiral  arrays  redirects 
a  larger  fraction  of  the  incident  radiation  into  the  absorbing  Si  substrate.  Moreover,  this 
effect  increases  the  optical  path  of  photons  in  the  photodetector,  as  well  as  enhancing  the 
coupling  to  LSPs  in  the  array  plane.  In  Fig.  10(c)  we  show  the  measured  photocurrent 
spectra  for  the  best  performing  GA  spiral  photodetector  device.  The  reference  line  shows 
the  photocurrent  spectra  measured  on  the  unpatterned  devices  in  the  nearest  reference 
cells,  respectively.  A  maximum  photocurrent  enhancement  of  approximately  a  factor  of  3 
is  measured  for  the  GA  spiral  at  950nm.  Such  an  increase  in  photocurrent  results  from  the 
interplay  of  plasmonic -photonic  coupling  in  GA  spirals,  which  contribute  to  the  overall 
enhancement  as  follows:  by  providing  more  efficient  coupling  of  incident  radiation  into 
the  thin  Si  layer  due  to  photonic  wide-angle  scattering;  by  enhancing  the  intensity  of  the 
near-fields  around  the  nanoparticles  at  the  Si  interface  owing  to  better  LSP  coupling  in 
the  plane  of  the  array  compared  to  periodic  structures.  A  detailed  comparison  with 
optimized  periodic  grating  structures  in  discussed  in  Ref.  [22]. 

In  Fig.  10(d)  we  show  the  integrated  photocurrent  enhancement  ratio,  calculated  by  the 
ratio  of  the  integrated  photocurrent  spectrum  of  the  device  with  and  without  plasmonic 
arrays.  The  ratios  falling  below  the  dotted  line  indicate  devices  with  overall  reduced 
performance  when  compared  against  their  neighboring  empty  reference  cells.  We  see  in 
Fig.  10(d)  that  both  GA  spiral  and  periodic  arrays  exhibit  an  optimization  trend  with 
respect  to  the  interparticle  spacing,  yielding  maximum  integrated  enhancements  of  8% 
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and  31%  over  the  reference  cells,  respectively.  We  observe  finally  that  the  experimentally 
measured  31%  integrated  enhancement  of  GA  spiral  arrays  has  been  demonstrated 
covering  only  25%  of  the  active  photodetector  area  by  GA  arrays  (Fig  10c,  inset). 

The  significant  photocurrent  enhancement  values  demonstrated  for  devices  coupled  to 
GA  spirals  highlight  the  potential  for  even  greater  performances  in  the  case  of  complete 
coverage  of  the  device  area. 

Spatial-spectral  detection  with  colorimetric  fingerprints 

In  this  project  we  have  developed  a  novel  concept  for  optical  biosensing  based  on  the 
engineering  of  colorimetric  fingerprints  in  aperiodically  nanopattemed  surfaces  [4,5]. 
Light  scattered  from  these  surfaces  is  localized  in  different  areas  of  the  devices  and  gives 
rise  to  complex  color  structures  that  specifically  fingerprint  the  aperiodic  surface. 

In  current  biosensing  technology,  two-dimensional  optical  gratings  provide  a  well- 
established  platform  for  biochemical  colorimetric  detection  and  enabled  label-free 
sensing  of  various  molecular  analytes  and  protein  dynamics.  Periodic  grating  biosensors 
provide  a  distinct  change  in  the  intensity  of  diffracted  light  or  in  the  frequency  spectrum 
of  optical  resonances  in  response  to  variations  of  refractive  index.  The  physical 
mechanism  at  the  basis  of  the  colorimetric  response  is  the  well-known  phenomenon  of 
Bragg  scattering.  While  this  process  provides  frequency  selective  responses  that  are 
useful  for  colorimetric  detection,  the  ability  of  light  waves  to  interact  with  adsorbed  or 
chemically  bound  analytes  present  on  the  surface  of  the  sensors  is  intrinsically  limited.  In 
fact,  Bragg  scattering  is  a  first-order  process  in  surface  scattering  perturbation  theory,  and 
scattered  photons  easily  escape  from  a  periodic  surface  without  prolonged  interaction 
with  the  sensing  layer.  On  the  other  hand,  engineering  multiple  light  scattering  in  planar 
aperiodic  structures  provides  novel  opportunities  for  bio-chemical  sensing  applications 
because  it  enables  much  longer  photon  dwelling  times  with  the  sensing  areas.  In 
particular,  optical  sensing  platforms  can  be  boosted  by  developing  scattering  elements 
that  simultaneously  provide  high  sensitivity  to  the  environmental  changes  and  high 
spectral  resolution,  as  both  factors  contribute  to  the  improvement  on  the  sensor  detection 
limit.  Detector  sensitivity  is  conventionally  defined  as  the  magnitude  of  the  wavelength 
shift  induced  by  a  change  in  the  ambient  refractive  index  (measured  in  nm/RIU),  and  can 
be  improved  by  enhancing  light-matter  interactions  on  a  planar  substrate.  In  turn,  the 
resolution  in  measuring  wavelength  shifts  inversely  depends  on  the  linewidth  of  the 
resonant  mode  supported  by  the  structure.  We  have  demonstrated  that  aperiodic  photonic 
structures  provide  the  necessary  balance  between  the  resonant  character  of  their  resonant 
modes  and  the  spatial  distribution  of  large  field  intensity  over  extended  sensing  areas, 
resulting  in  largely  improved  sensitivity  over  periodic  grating  sensors  and  or  photonic 
crystals  cavities,  which  are  limited  by  the  small  overlap  of  the  analyte  with  the  localized 
field.  Building  on  these  results,  we  have  introduced  within  this  prpoject  a  novel  approach 
to  label-free  optical  biosensing  based  on  micro-spectroscopy  and  spatial  correlation 
analysis  of  structural  color  patterns  excited  by  white  light  illumination  in  DANS  metal- 
dielectric  surfaces  [4].  In  contrast  to  traditional  photonic  gratings  or  photonic  crystals 
sensors  (which  efficiently  trap  light  within  small-volumes),  aperiodic  scattering  sensors 
sustain  distinctive  resonances  localized  over  larger  surface  areas. 
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Fig.  11.  (a,b)  Colorimetric  fingerprint  of  as  a2  aperiodic  spiral  array  of  A1  nanoparticles  measured  by  dark- 
field  microscopy  and  an  SEM  picture  (b)  showing  a  detail  of  the  structure.  SEM  pictures  of  (a)  Fibonacci 
array,  (d)  Rudin-Shapiro  array,  and  (e)  Gaussian  prime  array  of  Au  nano -cylinders  with  radius  r=100nm, 
height  h=30nm,  and  minimum  interparticle  separation  d=25nm.  The  color  figures  in  (c-e)  are  the 
corresponding  “colorimetric  fingerprints”  measured  by  dark-field  microscopy  under  white  light 
illumination. 

These  structural  resonances,  known  as  critical  modes,  are  described  by  highly  fluctuating 
field  profiles  described  by  multi-fractal  analysis  and  possess  a  dense  spectrum,  resulting 
in  efficient  photon  trapping  through  higher-order  multiple  scattering  processes  on  the 
surface  of  the  devices.  These  distinctive  scattering  behavior  strongly  enhance  the 
sensitivity  of  DANS  surfaces  to  small  refractive  index  variations.  The  complex  spatial 
patterns  of  critical  modes  offer  the  potential  to  engineer  structural  color  sensing  using 
spatially  localized  patterns  at  multiple  wavelengths.  In  fact,  when  illuminated  by  white 
light,  DANS  surfaces  feature  highly  structured  multi-color  scattering  patterns  which  are  a 
sensitive  fingerprints  of  the  aperiodic  geometries,  known  as  colorimetric  fingerprints. 
When  light  is  incident  on  aperiodic  surfaces  with  engineered  geometry  and  particle 
separation,  distinct  optical  frequencies  (i.e.,  colors)  localize  in  different  areas  of  the 
device  and  can  readily  be  observed  in  the  object-plane  or  in  the  far- zone,  as  demonstrated 
by  the  color  patterns  in  Fig.  11.  Adding  a  thin  layer  of  analyte  on  top  of  aperiodic  arrays 
shifts  the  resonance  wavelengths  of  the  resonances  supported  by  the  aperiodic  surface 
leading  to  a  spatial  rearrangement  of  the  localized  field  intensity  across  the  area.  Due  to 
the  critical  nature  of  the  resonant  states,  which  are  spatially  extended  though  highly 
fluctuating  in  space,  a  global  change  in  the  colorimetric  pattern  of  scattered  radiation  can 
readily  be  detected. 
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Fig.  12.  (a)  Coating  of  different  thicknesses  of  silk  protein  monolayers  were  characterized  by  Atomic  Force 
Microscopy  (insert)  and  the  colorimetric  responses  of  the  associated  arrays  were  measured  spectrally. 
Coating  of  different  thicknesses  of  silk  protein  monolayers  were  characterized  by  Atomic  Force 
Microscopy  (insert)  and  the  colorimetric  responses  of  the  associated  arrays  were  measured  spectrally,  (b) 
The  sensitivity  of  the  arrays  is  quantified  by  the  spectral  shift  of  the  scattered  radiation  peaks  (PWS)  per 
thickness  variation  of  the  protein  layer,  (d)  1-D  ACF  profiles  extracted  from  2-D  normalized 
autocorrelation  function  along  the  x-axis  of  the  middle  of  the  images,  (c)  The  changes  of  patterns  due  to 
different  thicknesses  of  silk  protein  monolayers  are  quantified  by  the  normalized  ACF  variances. 

The  colorimetric  fingerprints  of  nano-patterned  deterministic  aperiodic  surfaces  occur 
over  a  broad  spectral-angular  range  and  can  be  captured  using  conventional  dark-field 
microscopy,  as  demonstrated  in  Fig.  11.  Asa  result,  engineered  colorimetric  fingerprints 
of  aperiodic  surfaces  can  be  used  as  a  transduction  mechanism  in  a  novel  type  of  highly 
sensitive  label-free  multiplexed  sensor.  In  such  sensing  device,  both  the  peak  wavelength 
shifts  (Fig.  12a, b)  of  scattered  radiation  as  well  as  the  environment-dependent  spatial 
structure  of  the  colorimetric  fingerprints  (Fig.  12  c-d)  can  be  utilized  to  detect  the 
presence  of  nanoscale  protein  layers  on  the  surface  of  DANS.  The  proposed  approach  is 
intrinsically  more  sensitive  to  local  refractive  index  modifications  compared  to  traditional 
ones  due  to  the  enhancement  of  small  phase  variations,  which  results  from  the  multiple 
scattering  regime. 

The  spatial  modifications  of  the  structural  color  fingerprints  induced  by  small  refractive 
index  variations  can  be  quantified  using  standard  image  autocorrelation  analysis 
performed  directly  on  the  scattered  radiation.  This  can  be  accomplished  in  two 
modalities:  (a)  using  polychromatic  or  (b)  single  wavelength  radiation,  as  shown  in  Fig. 
12  in  the  case  of  silk  protein  monolayers  deposited  on  engineered  DANS  surfaces.  To 
construct  the  image  autocorrelation  function,  the  value  of  the  field  intensity  at  point  (x,  y) 
in  the  array  plane  is  compared  with  that  at  another  point  (x',  y')  and  mapped  as  a  function 
of  the  distance  between  the  two  points.  ID  profiles  through  the  2D  autocorrelation 
function  of  the  intensity  images  in  Fig.  12(d).  The  initial  decay  in  the  ACF  reflects  local 
short-range  correlations  in  the  aperiodic  structure,  and  any  long-range  periodicities  in  the 
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intensity  pattern  give  rise  to  periodic  oscillations  in  the  ACF.  To  quantify  the  overall 
change  in  the  colorimetric  pattern  induced  by  the  presence  of  thin  molecule  layers,  we 
calculate  the  variance  of  the  fluctuations  in  the  intensity  distribution  function,  which  can 
be  found  as  the  value  of  the  properly  normalized  discrete  ACF  in  the  limit  of  zero  lateral 
displacements.  The  proposed  approach  recently  led  to  the  demonstration  of  femto-molar 
detection  of  silk  proteins  adsorbed  on  DANS. 

The  refractive  index  modifications  induced  by  the  analytes  (silk  protein)  can  be  detected 
by  frequency  shifts  and  a  global  structural  color  modification.  Combining  Electron  Beam 
Lithography  (EBL),  dark-field  scattering  micro-spectroscopy,  autocorrelation  analysis 
and  rigorous  multiple  scattering  calculations  based  on  the  Generalized  Mie  Theory 
(GMT).  In  this  project  we  have  engineered  aperiodic  arrays  of  metal  nano-particles  on 
quartz  substrates,  and  showed  that  the  information  encoded  in  both  the  spectral  and 
spatial  distribution  of  structural  resonances  can  be  simultaneously  utilized  for  sensitive 
bio-detection  on  the  nanoscale  [4,5,28].  The  potential  of  the  proposed  approach  for  rapid, 
label-free  detection  and  recognition  of  biomolecular  analytes  in  the  visible  spectral  range 
was  experimentally  demonstrated  by  the  distinct  variation  in  the  spectral  and  spatial 
colorimetric  fingerprints  in  response  to  monolayer  increments  of  protein  layers 
sequentially  deposited  on  the  surface  of  aperiodic  arrays  of  nanoparticles. 

Recently,  we  also  demonstrated  successful  integration  of  aperiodic  arrays  of  metal 
nanoparticles  with  microfluidics  technology  for  optical  sensing  using  the  spectral- 
colorimetric  responses  of  the  nanostructured  arrays  to  refractive  index  variations  (S.Y. 
Lee  et  al.,  in  review).  In  particular,  they  fabricated  different  aperiodic  arrays  of  Au 
nanoparticles  with  varying  interparticle  separations  and  Fourier  spectral  properties  using 
Electron  Beam  Lithography  (EBL)  and  integrated  them  with  polydimethylsiloxane 
(PDMS)  microfluidics  structures  by  soft-lithographic  micro-imprint  techniques. 

The  fabrication  process  flow  that  enabled  the  optofluidic  integration  of  DANS  is 
schematically  illustrated  in  Fig.  13  (a).  In  our  integration  process,  a  28pm-thick  SU8 
negative  photoresist,  with  an  additional  lpm-thick  SU8  filling  groove  for  facilitating 
fluid  filling  into  microchannels,  was  first  spin-coated  on  a  cleaned  silicon  wafer  and 
microfluidic  structures  were  defined  on  the  substrate  through  standard  photolithography. 
A  thin  layer  of  PDMS  (10:1  base:  curing  agent)  was  subsequently  cured  on  top  of  the 
SU8  master  at  70°C  for  2  hours.  Half-through  inlet  and  outlet  holes  were  bored  on  the 
PDMS  mold  at  both  ends  of  the  transferred  microchannel  patterns  and  the  PDMS  mold 
was  oxygen-plasma  treated  along  with  the  colorimetric  DANS  sensor  on  SiOo  substrate  to 
form  an  optofluidic  DANS  device. 

The  oxygen  plasma  process  permanently  bonded  the  device  together  and  resulted  in 
hydrophilic  surfaces,  which  easily  prime  the  microchannels  with  insensible  amount  of 
pressure.  Epidermic  syringe  needles  were  then  inserted  from  the  sides  of  the  PDMS  mold 
into  the  holes  as  fluid  access  ports  and  secured  using  silicone  sealant,  as  shown  in  Fig.  13 
(b,c).  Figure  13(b)  clearly  illustrates  optofluidic  DANS  devices  under  a  microscope,  in 
which  various  types  of  arrays  were  positioned  inside  the  microfluidic  flow  cell  for 
colorimetric  detection. 
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Fig.  13.  (a)  Fabrication  process  flow  of  PDMS  micro-channel  and  microfluidic  integration  with 
colorimetric  DANS  sensing  structures,  (b)  optical  microscope  image  of  the  microfluidic  flow  cell  with 
different  2D  aperiodic  arrays.  The  scale  bar  equals  100pm.  (c)  Digital  image  of  a  packaged  opto  fluidic 
DANS  device,  (d)  Dark-field  microscope  setup  utilized  for  sensing  characterization  of  fabricated  devices, 
(e)  Plasmon-enhanced  linear  scattering  shift  of  microfluidics  integrated  DANS  substrates,  (g)  Variation  of 
autocorrelation  functions  of  the  scattered  radiation  from  microfluidics  integrated  DANS  sensors  immersed 
in  a  liquid  with  varying  refractive  index,  (h)  Measured  dependence  of  the  autocorrelation  functions  of  the 
scattered  radiation  in  a  microfluidics  channel  as  a  function  of  the  liquid  refractive  index. 

Spectral  shifts  of  scattering  spectra  and  the  distinctive  modifications  of  structural  color 
patterns  induced  by  refractive  index  variations  were  simultaneously  measured  inside 
microfluidic  flow  cells  by  dark-field  spectroscopy  and  image  correlation  analysis  in  the 
visible  spectral  range,  as  demonstrated  in  Figs.  13  (e-h).  The  main  outcome  of  this  study 
is  that  aperiodic  plasmonic  Galois  arrays  result  in  large  fluctuations  of  colorimetric 
fingerprints  simultaneously  accompanied  by  large  wavelength  shifts  of  their  scattering 
spectra,  providing  both  high  spectral  and  colorimetric  sensitivities  within  a  microfluidic 
device.  The  integration  of  plasmonic  aperiodic  arrays  with  microfluidics  devices  provides 
a  novel  approach  with  multiplexed  spatial-spectral  responses  for  opto-fluidics  platforms 
and  lab-on-a-chip  optical  biosensing. 

Engineering  optical  angular  momentum  of  light  with  aperiodicity 

During  this  project  we  have  discovered  that  Vogel  spiral  nanoparticle  arrays,  when 
illuminated  by  optical  beams,  give  rise  to  scattered  radiation  carrying  OAM  [8].  Dal 
Negro  et  al.  [27]  developed  an  analytical  model  that  captured  in  closed  form  solution  the 
Fourier  and  diffraction  properties  of  arbitrary  Vogel  spiral  arrays  and  demonstrated  that 
light  diffraction  by  Vogel  spirals  encode  in  the  far-field  OAM  values  arranged  in 
aperiodic  sequences  associated  to  number-theoretic  properties  of  the  spirals.  More 
precisely,  the  OAM  values  diffracted  in  the  far-field  are  determined  by  the  rational 
approximations,  given  by  the  continued  fraction  expansion,  of  the  irrational  divergence 
angles  of  Vogel  spirals.  In  particular,  wave  diffraction  by  GA  arrays  encodes  a  Fibonacci 
sequence  of  OAM  values  in  the  Fraunhofer  far  field-region  [27]. 
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Fig.  14.  (a,b)  (a,d)  SEMs  of  GA  and  (a.  spirals  respectively  composed  of  500nm  diameter  gold  nanoparticles 
on  fused  silica.  (c,d)  Analytically  calculated  far  field  intensity  from  the  spirals  in  (a,b)  at  633nm.  (e) 
Optical  set  up  used  for  the  complex  amplitude  retrieval  (f,h)  Analytically  calculated  phase  portrait  of  far 
field  radiation  from  GA  and  p  spirals  respectively.  (g,i)  Experimentally  measured  phase  portrait  from  GA 
and  p  spirals,  respectively.  Shaded  areas  indicate  regions  where  the  field  intensity  is  too  low  to  measure  the 
phase  accurately.  These  regions  are  ignored  in  the  FHD. 

Controlled  generation  and  manipulation  of  OAM  states  with  large  values  of  azimuthal 
numbers  has  been  recently  demonstrated  experimentally  in  this  project  by  Lawrence  et  al. 
[31]  using  various  types  of  Vogel  arrays  of  metal  nanoparticles.  The  OAM  content  of  a 
diffracted  optical  wave  at  633nm  was  analyzed  using  phase  stepped  interferometric 
measurements  to  recover  the  complex  optical  field.  Fourier-Hankel  modal  decomposition 
of  scattered  radiation  was  performed  to  demonstrate  the  generation  of  OAM  sequences 
from  Vogel  spirals.  The  samples  consisted  in  arrays  of  Au  nanoparticles  fabricated  on 
fused  silica  using  electron  beam  lithography,  metal  evaporation  and  lift-off  processing 
[31].  The  SEMs  of  fabricated  arrays  are  shown  in  Figs.  14(a,b).  The  corresponding  far- 
field  diffraction  patterns  are  shown  in  Figs.  14  (c,d).  In  Fig.  14(e)  we  show  the  optical 
set-up  used  to  measure  the  complex  far  field.  A  HeNe  laser  was  weakly  focused  to  a  spot 
size  of  50pm,  from  the  rear  of  the  sample  and  the  scattered  light  from  the  array  was 
collected  by  a  50X  (NA=.75)  objective.  Two  additional  lenses  were  added  to  image  the 
far  field  at  the  plane  of  the  CCD[31]. 
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In  order  to  measure  the  phase  of  the  far  field  radiation  a  reference  beam  was  reflected 
from  the  piezo  stage  mounted  mirrors,  expanded  and  directed  to  the  CCD.  Multiple 
interference  patterns  were  collected  by  increasing  the  piezo  bias  voltage,  scanning  the 
phase  of  the  reference  beam.  A  phase  retrieval  algorithm  was  finally  used  to  recover  the 
phase  of  the  scattered  light  relative  to  the  one  of  the  reference  beam.  The  analytically 
calculated  phase  of  the  far  field  radiation  from  GA  and  p  spirals  is  shown  in  Figs.  14(f,h) 
respectively.  The  measured  phase  from  a  GA  spiral  is  shown  in  Fig.  14(g)  and  the  phase 
from  a  p  spiral  is  shown  in  Fig.  22(i). 

Modal  decomposition  was  used  to  analyze  a  superposition  of  OAM  carrying  modes  in  the 
far  field  pattern  and  determine  their  relative  contribution.  Decomposition  into  a  basis  set 
with  azimuthal  dependence,  is  accomplished  through  Fourier-Hankel  decomposition 
(FHD)  according  to  [15,18,25,27,31]: 
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where  Jm  is  the  m-th  Bessel  function.  In  this  decomposition  the  m-th  order  function 
carries  OAM  with  azimuthal  number  m,  this  can  accommodate  positive  and  negative 
integer  values  for  m.  Since  we  are  primarily  concerned  with  the  azimuthal  component  of 

-p  f  yyi  \ 

the  field  we  can  sum  J  v  ’  r  over  radial  wavevectors  kr.  In  Figs.  15(a,b)  the  summed 
FHD  of  the  calculated  far  field  radiation  from  GA  and  p  spirals  are  shown,  respectively. 
Both  show  peaks  at  azimuthal  numbers  corresponding  to  the  denominators  of  the  rational 
approximations  (i.e.,  convergents)  of  the  irrational  divergence  angles  used  to  generate  the 
spirals,  as  amply  discussed  in  Ref.  [27].  For  the  GA  spiral  this  results  in  encoding  the 
Fibonacci  sequence  in  the  OAM  spectrum.  The  p  spiral  has  significantly  fewer 
convergents  in  our  region  of  interest,  however  some  OAM  peaks  are  observed  at  linear 
combinations  of  convergents.  These  peaks  can  be  reduced  with  aperturing  of  the  beam. 
The  values  of  OAM  peaks  can  be  explained  based  on  the  analytical  far  field  radiation 
pattern  theory  developed  by  Dal  Negro  [27].  Here  it  suffices  to  remind  that  the  FHD  of 
the  Vogel  far  field  radiation  pattern,  is  given  by  [27]: 

f(m,kr)  =  f,Am(kry™ 

n-\ 


A  (k  ) 

where,  N  is  the  number  of  particle  in  the  array  and  m  v  '  gives  kr  dependent  coefficients 
which  we  will  ignore  as  we  are  primarily  concerned  with  the  azimuthal  component. 
When  ma  is  an  integer,  the  N  waves  in  the  equation  above  will  be  in  phase  and  produce 
an  OAM  peak  at  azimuthal  number  m.  If  a  is  an  irrational  number,  this  condition  will 
never  be  satisfied.  However,  an  irrational  can  be  approximated  by  rational  fractions 
(convergents)  which  make  the  product  ma  approximately  an  integer  when  m  equals  the 
different  denominators  of  the  convergents  [27]. 
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Theory 


Fig.  15.  Fourier-Hankel  decomposition,  summed  over  radial  wavevector,  for  GA  spirals  (a,c)  and  p  spirals 
(b,d).  In  (a,b)  FHD  is  obtained  using  analytically  calculated  far  field,  while  in  (c,d)  FHD  is  obtained  using 
experimentally  measured  far  field. 

This  result  was  demonstrated  by  examining  the  OAM  content  of  the  experimentally 
measured  scattered  radiation.  The  reconstructed  complex  field  (Figs.  14(g,i)),  is 
decomposed  into  azimuthal  and  radial  components.  The  shaded  areas  in  Figs.  14(f-i) 
indicate  the  limits  (angular  range)  where  the  measurement  of  the  scattered  light  is 
accurate  due  to  the  limitations  of  the  experimental  setup.  Figure  15  directly  compares  the 
predicted  OAM  spectra  (i.e.,  radially  summed  FHD)  and  the  experimentally  measured 
ones.  In  Fig.  15(c)  experimental  results  are  shown  for  the  GA  spiral  and  peaks  in  the 
transform  are  observed  at  numbers  corresponding  to  the  Fibonacci  sequence.  In  Fig. 
15(d)  peaks  are  observed  matching  the  theoretically  predicted  values  for  the  p  spiral.  A 
full  list  of  analytically  calculated  and  the  experimentally  measured  FHD  peak  positions 
can  be  found  in  Ref  [31]. 

Finally,  we  have  demonstrated  direct  excitation  of  optical  modes  carrying  OAM  values  in 
light-emitting  dielectric  GA  spiral  structures  consisting  of  Er-doped  SiN  nanopillars  []. 
The  successful  demonstration  of  the  superposition  of  many  OAM  modes  in  engineered 
Vogel  spiral  arrays  of  Au  nanoparticles  provides  novel  exciting  opportunities  for 
radiation  engineering  that  are  of  interest  to  secure  optical  communication,  classical,  and 
quantum  cryptography. 
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Deterministic  aperiodic  lasers 

In  this  project  we  finally  demonstrated,  in  collaboration  with  Yale  University,  the  first 
laser  action  from  a  multiply  scattering  pseudo-random  medium  [7].  This  is  a 
breakthrough  results  in  the  field  of  complex  lasers  since  we  have  experimentally 
demonstrated  that  deterministic  aperiodic  lasers  behave  similarly  to  random  lasers  but  the 
frequency  positions  of  their  many  lasing  modes  are  completely  predictable. 


Fig.  16.  Cover  of  Applied  Physics  Letters  “ Demonstration  of  laser  action  in  a  pseudorandom  medium ”, 
97,223101  (2010) 

We  have  demonstrated  experimentally  that  the  localized  lasing  modes  in  the  Rudin- 
Shapiro  arrays  of  air  nano-holes  in  GaAs  membranes  occur  at  reproducible  spatial 
locations  and  their  frequencies  are  only  slightly  affected  by  the  structural  fluctuations  in 
different  samples.  Numerical  study  on  the  resonances  of  the  passive  systems  and  optical 
imaging  of  lasing  modes  enabled  us  to  interpret  the  observed  lasing  behavior  in  terms  of 
distinctive  localized  resonances  in  the  two-dimensional  Rudin-Shapiro  (RS)  structures. 
Although  the  relatively  narrow  gain  spectrum  of  InAs  quantum  wells  does  not  allow 
broad-band  lasing  in  the  current  samples,  our  numerical  simulations  confirmed  that  the 
RS  structures  support  many  localized  modes  of  distinct  frequencies.  We  expect  lasing 
over  a  wide  frequency  range  once  switching  to  broad-band  gain  media  such  as 
semiconductor  quantum  dots.  The  deterministic  aperiodic  Rudin-Shapiro  systems 
provides  an  alternative  approach  from  random  media  and  photonic  crystals  for  the 
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engineering  of  multi-frequency  coherent  light  sources  and  complex  cavities  amenable  to 
predictive  theories  and  technology  integration. 
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Fig.  17:  (a)  Emission  spectra  taken  at  the  incident  pump  power  P=  3,  4.1,  6,  9,  15,  20  pW.  (b)  Measured 
intensity  (black  circle)  and  spectral  width  (blue  square)  of  one  emission  peak  at  wavelength  1=  978  nm  vs. 
P  showing  the  onset  of  lasing  at  ~  3  pW.  The  red  lines  are  linear  fitting  of  peak  intensity  below  and  above 
the  lasing  threshold  of  the  device. 

In  addition,  we  demonstrated  lasing  action  in  quasi-2D  Thue-Morse  structures  with 
optimal  degree  of  aperiodicity  fabricated  in  a  GaAs  membrane  structure  [10].  By 
changing  the  relative  size  of  two  scatterers  that  correspond  to  the  building  blocks  A  and  B, 
we  gradually  vary  the  degree  of  aperiodicity  and  investigate  its  effect  on  lasing  (Fig.  17). 
We  discovered  that  there  exists  an  optimal  degree  of  aperiodicity  where  the  quality  factor 
reaches  a  global  maximum  and  lasing  becomes  the  strongest.  This  is  attributed  to  an 
enhancement  of  horizontal  confinement  of  light  in  a  finite-size  pattern  by  structural 
aperiodicity.  However,  the  continuous  background  in  the  spatial  Fourier  spectrum  of  the 
Thue-Morse  structure  facilitates  vertical  leakage  of  light  out  of  the  membrane.  At  various 
degrees  of  aperiodicity,  different  types  of  modes  acquire  the  highest  quality  factors  and 
may  be  selected  for  lasing.  This  work  opens  a  way  of  controlling  lasing  characteristic  via 
tuning  structural  aperiodicity. 


(a)  (b) 
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